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Abstract. We review 1 a few off-the-beaten-track ideas in cosmology. They 
solve a variety of fundamental problems; also they are fun. We start with 
a description of non-singular dilaton cosmology. In these scenarios gravity 
is modified so that the Universe does not have a singular birth. We then 
present a variety of ideas mixing string theory and cosmology. These solve 
the cosmological problems usually solved by inflation, and furthermore shed 
light upon the issue of the number of dimensions of our Universe. We finally 
review several aspects of the varying speed of light theory. We show how the 
horizon, flatness, and cosmological constant problems may be solved in this 
scenario. We finally present a possible experimental test for a realization of 
this theory: a test in which the Supernovae results are to be combined with 
recent evidence for redshift dependence in the fine structure constant. 



1. Introduction 

In spite of their unprecedented success at providing causal theories for the 
origin of structure, our current models of the very early Universe, in partic- 
ular models of inflation and cosmic defect theories, leave several important 
issues unresolved and face crucial problems (see [1] for a more detailed dis- 
cussion). The purpose of this chapter is to present some imaginative and 
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speculative ideas for early Universe cosmology which go beyond our present 
models and which address some of the puzzles unresolved by inflationary 
cosmology. 

Let us recall four of the key conceptual problems of inflationary cosmol- 
ogy The first is the fluctuation problem: without finely adjusting param- 
eters in the particle physics model, the amplitude of the resulting cosmic 
microwave background (CMB) anisotropies is generically predicted to be 
far in excess of the measured values [2]. What is missing is an in-built 
suppression mechanism for scalar metric fluctuations. 

The second key problem of potential-driven inflation is the super-Planck 
scale physics problem: unless the period of inflation is very short, comoving 
length scales corresponding to today's large scale structure had a physical 
wavelength much smaller than the Planck length in the early stages of 
inflation. In the absence of an implementation of inflation in the context 
of a unified theory of Planck-scale physics, it is not justified to blindly 
extrapolate the predictions of simple scalar field toy models to the required 
small length scales. 

Inflationary cosmology leaves the initial singularity problem of standard 
cosmology unresolved [3], and in this sense remains an incomplete theory 
of the early Universe without a consistent initial value formalation. This is 
a further problem of inflation. 

The Achilles heel of potential-driven inflationary cosmology is without 
doubt the cosmological constant problem. The mechanism of inflation relies 
on a transient cosmological constant generating the required exponential 
expansion of space, but can we trust this mechanism in the absence of an 
understanding of what sets the present value of the cosmological constant 
to (almost) zero? 

In the following, we summarize some imaginative but very speculative 
ideas which address one or more of the above problems. At the present 
time, none of the ideas is well enough developed to be able to solve all of 
the problems and none forms a complete well-developed theory. Our hope, 
however, is that the ultimate theory of the early Universe will contain cer- 
tain aspects of the mechanisms discussed below. As bi-products, these ideas 
point to the exciting possibility that the unified theory of all forces may 
explain the dimensionality of space [4], and yield an alternative to infla- 
tion for solving the flatness and horizon problems of standard cosmology 
without the need to invoke inflation [5, 6]. 

The outline of this chapter is as follows. In the next section, we show 
that adding a special class of higher derivative gravity terms to the usual 
low energy effective action for matter and gravity can lead to a nonsingu- 
lar Universe. Applied to dilaton gravity, this construction may lead to a 
graceful exit in "pre-big-bang" cosmology. 
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Although higher derivative gravity corrections to the Einstein action 
are a generic feature of any lower energy limit of a unified thoery, the 
specific terms used in the nonsingular Universe construction of Section 2 are 
rather ad hoc and not at the present time motivated by any fundamental 
microscopic theory. In Section 3 we study a few cosmological aspects of 
the one well-motivated candidate for a unified theory of all forces, string 
theory. In particular, we explore the role of target space duality in resolving 
cosmological singularities and in determining a preferential dimensionality 
of space. 

Perhaps the craziest attempt to bypass inflation is the varying speed 
of light cosmology [5, 6]. In its current formulation the theory allows for 
violations of Lorentz invariance. These are shown to lead to a solution of 
the problems usually solved by inflation, even if the matter content of the 
Universe satisfies the strong energy condition. We outline some aspects of 
this theory in Section 4, and in Section 5 we show how solutions to the 
cosmological problems emerge in these theories. In Section 6 we then show 
how such a theory may be confronted with experiment. 

In a concluding section we digress on the future of imaginative cosmol- 
ogy (as opposed to mainstream, "boring", cosmology). 

2. Nonsingular Dilaton Cosmology 

2.1. MOTIVATION 

The initial singularity is one of the outstanding problems of current cosmo- 
logical models. In standard big bang cosmology, the existence of the initial 
singularity is an inevitable consequence of the Penrose-Hawking theorems 
[7, 8]. While potential-driven inflationary models such as chaotic inflation 
[9] resolve many of the problems faced by conventional cosmology, initial 
singularities are still generic, even when stochastic effects are included [3]. 

It is generally hoped that string theory will lead to a resolution of the 
singularity problem. A key feature of string theory is that the low energy 
degrees of freedom include besides the graviton the dilaton and an antisym- 
metric tensor field. Gasperini and Veneziano initiated a program known as 
"pre-big-bang cosmology" [10] aimed at studying the cosmological conse- 
quences of the low energy effective action of string theory. To lowest order, 
and neglecting the antisymmetric tensor field, the action is given by 

where <f> is the dilaton, k 2 = 8itG = Sirm^ 2 , with G being the (4 dimen- 
sional) gravitational coupling, and m p i the Planck mass. The field equations 
of pre-big-bang cosmology exhibit a new symmetry, scale factor duality, 
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which (in the Einstein frame) maps an expanding Friedmann-Robertson- 
Walker (FRW) cosmology to a dilaton-dominated contracting inflationary 
phase. This raises the hope that it is possible to realize a nonsingular cos- 
mology in which the Universe starts out in a cold dilaton-dominated con- 
tracting phase, goes through a bounce and then emerges as an expanding 
FRW Universe (see Ref. [11] for a recent review of pre-big-bang cosmology 
Unfortunately it has been shown that the two branches of pre-big-bang cos- 
mology cannot be smoothly connected within the above tree-level action 
[12, 13, 14, 15]. The contracting dilaton-dominated branch has a future 
singularity, whereas the expanding branch emerges from a past singularity. 

A natural approach to resolving the singularity problem of general rel- 
ativity is to consider an effective theory of gravity which contains higher 
order terms, in addition to the Ricci scalar of the Einstein action. This 
approach is well motivated, since we expect that any effective action for 
classical gravity obtained from string theory, quantum gravity or by inte- 
grating out matter fields, will contain higher derivative terms. Thus, it is 
quite natural to consider higher derivative effective gravity theories when 
studying the properties of space-time at large curvatures. 

Most higher derivative gravity theories have much worse singularity 
problems than Einstein's theory. However, it is not unreasonable to expect 
that in the fundamental theory of nature, be it string theory or some other 
theory, the curvature of space-time is limited. In Ref. [16] the hypothe- 
sis was made that when the limiting curvature is reached, the geometry 
must approach that of a maximally symmetric space-time, namely de Sit- 
ter space. The question now becomes whether it is possible to find a class 
of higher derivative effective actions for gravity which have the property 
that at large curvatures the solutions approach de Sitter space. A nonsin- 
gular Universe construction which achieves this goal was proposed in Refs. 
[17, 18]. It is based on adding to the Einstein action a particular com- 
bination of quadratic invariants of the Riemann tensor chosen such that 
the invariant vanishes only in de Sitter space-times. This invariant is cou- 
pled to the Einstein action via a Lagrange multiplier field in a way that 
the Lagrange multiplier constraint equation forces the invariant to zero at 
high curvatures. Thus, the metric becomes de Sitter and hence explicitly 
nonsingular. 

If successful, the above construction will have some very appealing con- 
sequences. Consider, for example, a collapsing spatially homogeneous Uni- 
verse. According to Einstein's theory, this Universe will collapse in a finite 
proper time to a final "big crunch" singularity (top left Penrose diagram 
of Figure 1). In the new theory however, the Universe will approach a de 
Sitter model as the curvature increases. If the Universe is closed, there will 
be a de Sitter bounce followed by re-expansion (bottom left Penrose di- 
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agram in Figure 1). Similarly, spherically symmetric vacuum solutions of 
the new equations of motion will presumably be nonsingular, i.e., black 
holes would have no singularities in their centers. The structure of a large 
black hole would be unchanged compared to what is predicted by Einstein's 
theory (top right, Figure 1) outside and even slightly inside the horizon, 
since all curvature invariants are small in those regions. However, for r — > 
(where r is the radial Schwarzschild coordinate), the solution changes and 
approaches a de Sitter solution (bottom right, Figure 1). This would have 
interesting consequences for the black hole information loss problem (in two 
dimensions, this construction has been successfully realized [19]). 

In this section, we review the nonsingular Universe construction of [17, 
18] and its applications to dilaton cosmology [20, 21]. 

2.2. AN ANALOGY 

To motivate the nonsingular Universe construction, we turn to a well known 
example of a theory with a limited physical quantity, namely point particle 
motion in the theory of special relativity. The transition from the New- 
tonian theory of point particle motion to the special relativistic theory 
transforms a theory with no bound on the velocity into one in which there 
is a limiting velocity, the speed of light c (in the following we use units in 
which h = c = 1). This transition can be obtained [17] by starting with the 
Newtonian action of a point particle with world line x(t): 



and adding a Lagrange multiplier ip which couples to x 2 , the quantity to 
be made finite, and which has a potential V(tp): 



From the Lagrange multiplier constraint equation 



it follows that x 2 is limited provided V((p) increases no faster than lin- 
early in ip for large \(p\. The small ip asymptotics of V((p) is determined by 
demanding that at low velocities the correct Newtonian limit results: 




(2) 




(3) 



V((p) ~ f 2 
V(<f) ~ ip 



as\(p\ — > , 
as \<p\ — ► oo . 



(5) 
(6) 
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Choosing the simple interpolating potential 

the Lagrange multiplier can be integrated out, resulting in the well-known 
action 

S„ cw = \ j dtVl ~ X 2 (8) 

for point particle motion in special relativity. 



2.3. CONSTRUCTION 

The procedure for obtaining a nonsingular Universe theory [17] is based on 
generalizing the above Lagrange multiplier construction to gravity. Start- 
ing from the Einstein action, one can introduce a Lagrange multiplier <pi 
coupled to the Ricci scalar R to obtain a theory with limited R: 

s = J d 4 xV=g(R + viR + Vi(vi)) , (9) 

where the potential Vi(<pi) satisfies the asymptotic conditions (5,6). 

However, this action is insufficient to obtain a nonsingular gravity the- 
ory. For example, singular solutions of the Einstein equations with R = 
are not affected at all. The minimal requirements for a nonsingular the- 
ory are that all curvature invariants remain bounded and the space-time 
manifold is geodesically complete. Implementing the limiting curvature hy- 
pothesis [16], these conditions can be reduced to more manageable ones. 
First, we choose one curvature invariant Iiig^v) (e.g. I\ = R in (9)) and de- 
mand that it be explicitely bounded by the construction of (9). In a second 
step, we demand that as h(g^ u ) approaches its limiting value, the metric 
g^v approach the de Sitter metric , a definite nonsingular metric with 
maximal symmetry. In this case, all curvature invariants are automatically 
bounded (they approach their de Sitter values), and the space-time can be 
extended to be geodesically complete. 

One can implement the second step of the above procedure by an- 
other Lagrange multiplier construction [17]. Consider a curvature invariant 
hig^u) with the property that 

h(g^) = o e> g^ = g^ s . (10) 

Next, introduce a second Lagrange multiplier field (f2 which couples to I2 
and choose a potential V 2 {^>2) which forces I2 to zero at large \<P2\- 



S = / d 4 x^g-[R + + Vifa) + 992/2 + V 2 (<P2)] , (11) 
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with asymptotic conditions (5,6) for V\_((pi) and conditions 



^2(^2) 



const as | ip 2 \ — > 00 



ip 2 as \<p 2 \ 



0, 



(12) 
(13) 



for V 2 (ip 2 ). The first constraint forces I 2 to zero, the second is required in 
order to obtain the correct low curvature limit. 

These general conditions are reasonable, but not sufficient in order to 
obtain a nonsingular theory. It must still be shown that all solutions arc 
well behaved, i.e., that they asymptotically reach the regions | <^>2 1 — > 00 of 
phase space (or that they can be controlled in some other way). This must 
be done for a specific realization of the above general construction. 

At the moment we are only able to find an invariant I 2 which singles out 
de Sitter space by demanding I 2 = provided we assume that the metric 
has special symmetries. The choice 



h 



R 2 + C 2 f/\ 



(14) 



singles out the de Sitter metric among all homogeneous and isotropic met- 
rics (in which case adding C 2 , the Weyl tensor square, is superfluous), all 
homogeneous and anisotropic metrics, and all radially symmetric metrics. 
As a specific example one can consider the action [17, 18] 



S = J d A x^g~ R + ^R- {<p 2 + 4^iH 2 1/2 + + V2&2 



with 



12 H, 



2 fi 



1 



ln(l + <pi) 



1 + <Pi V 1 + <Pi 

2 

1 -f ip 2 



(15) 

(16) 
(17) 



When restricted to the set of homogeneous and isotropic space-times, 
the equations of motion which follow from the action (15) are not hard 
to analyze [17, 18]. They consist of the two Lagrange multiplier constraint 
equations and the goo equation. The Lagrange multiplier constraint equa- 
tions can be used to eliminate the Hubble parameter and its derivative in fa- 
vor of functions of the Lagrange multiplier fields. The remaining one dynam- 
ical equation becomes an equation for the flow in the (ipi,ip 2 ) phase space. 
It can be shown that all solutions are either periodic about Minkowski 
space-time (<pi, <p 2 ) = (0, 0) or else asymptotically approach de Sitter space 
(\(p 2 \ -► 00). 
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One of the most interesting properties of this theory is asymptotic free- 
dom [18], i.e., the coupling between matter and gravity goes to zero at high 
curvatures. It is easy to add matter (e.g., dust, radiation or a scalar field) to 
the gravitational action in the standard way. One finds that in the asymp- 
totic de Sitter regions, the trajectories of the solutions projected onto the 
(^l) ¥>2) plane are unchanged by adding matter. This applies, for example, 
in a phase of de Sitter contraction when the matter energy density is in- 
creasing exponentially but does not affect the metric. The physical reason 
for asymptotic freedom is obvious: in the asymptotic regions of phase space, 
the space-time curvature approaches its maximal value and thus cannot be 
changed even by adding an arbitrarily high matter energy density. 

In retrospect, it is not necessary to explicitly bound an invariant I\ to 
obtain a nonsingular Universe. In fact, (p± can be eliminated from the ac- 
tion (15) without affecting the nonsingularity of the resulting cosmological 
solutions, as was explicitly demonstrated in [18]. Hence, in our application 
to dilaton cosmology we will only keep a single Lagrange multiplier field (f2 
which we will henceforth denote by ip. 

2.4. APPLICATION TO DILATON COSMOLOGY 

The starting point is the action (1) for dilaton gravity (written in the 
Einstein frame) to which we [20] add the higher derivative term given by 
I2, in analogy to what was done above in the absence of the dilaton. To be 
specific, we choose minimal coupling of I2 to the dilaton. 

Restricted to a homogeneous and isotropic metric of the form 

ds 2 =dt 2 -a(t) 2 ( L-dr 2 + rW), (19) 

1 — kr z 

where dQ, 2 is the metric on S 2 , the equations of motion resulting from (18) 
become 

= _ 3f ty +m - l(l x 2 + n 

H = -Vty)', 

X = -3H X , (20) 

where dots denote derivatives with respect to time, t, \ = <j> an d a prime 
(/) signifies the derivative with respect to ip. Since the initial goal is to 
construct a spatially flat, bouncing Universe, we have set the curvature 
constant k = 0. 
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Consider next the criteria which the potential V(ip) must satisfy. At 
small curvatures, the terms in the action (18) which depend on ip must be 
negligible compared to the usual terms of dilaton gravity. This is ensured 
by demanding 

Vty) ~ ip 2 M-0 (21) 

as the region of small \ip\ will correspond to the low curvature domain [18]. 
In order to implement the limiting curvature hypothesis, the invariant I2 
must tend to zero, and thus the metric g^ v will tend to a de Sitter metric 
at large curvatures, i.e. for \i\)\ — > 00. From the variational equation with 
respect to ip, it is obvious that this requires 

V(ip) -> const |V>|->oo. (22) 

Conditions (21) and (22) are the same as those (12,13) used in the pre- 
vious subsection, but they do not fully constrain the potential. In order to 
obtain a bouncing solution in the absence of spatial curvature it is neces- 
sary to add a third criterion: the equations must allow a configuration with 
H = and ip 7^ 0. From the equation of motion for ip in (20) it follows that 
V(ip) must become negative, assuming that it is positive for small \ip\. Let 
ipb denote the nontrivial zero of V(ip): 

%)=0. (23) 

In the absence of the dilaton, ipb will correspond to the value of ip at the 
bounce. In the presence of (p, the value of \ip\ at the bounce will depend 
on x an d will be larger than \ipb\. A simple potential which satisfies the 
conditions (21), (22) and (23) is [20] 



Vty) = t 16^_. (24) 

1 + ^V 4 



The conditions, (21) - (23), on the potential V(ip) discussed in the pre- 
vious section are necessary but not sufficient to obtain a nonsingular cos- 
mology. These conditions ensure that all solutions which approach large 
values of \ip\ are nonsingular, but the possibility of geodesic incompleteness 
for solutions which always remain within the small \ip\ region remains to be 
studied. If % = then the trajectories of the solutions in the phase plane 
(ip, H) can be studied quite easily analytically [20] and one can thereby 
explicitly demonstrate the absence of singularities. As a new feature (and 
as required in pre-big-bang cosmology), the theory admits spatially flat 
bouncing solutions. 

There are several special points and curves on the phase plane (ip,H). 
First, the point (ip, H) = (0, 0) corresponds to Minkowski space-time. The 
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potential V(tp) vanishes at this point, but it also vanishes at the points 

A = ±4. (25) 

Trajectories which pass through the phase plane points (ipj,, 0) correspond to 
spatially flat bouncing cosmologies. To demonstrate that the point (ip, H) = 
(4, 0) is in fact a bounce, one expands the ip equation of motion near H = 0, 
which yields 

Hip = —V . (26) 

When crossing the H = axis, the sign of ip changes. Contracting solutions 
with 2 < ip < 4 have ip > and approach the point (4, 0) in finite time 
since H is positive and does not tend to zero, and emerge for H > 
as expanding trajectories with decreasing curvature (since ip < 0). The 
trajectories in the phase plane are symmetric about the H = axis, except 
that the time arrows are reversed. 

By expanding the equations near the origin of the phase plane one can 
demonstrate the existence of critical lines, lines with 

located at 

ip c (H) = ±V6H (28) 

(for large values of \H\ the critical lines approach the lines ip = 2). Focusing 
on the contracting solutions, trajectories which lie above the critical line 
have H < and ip > 0. These trajectories thus are directed towards the 
line ip = 2 where H changes sign. Provided they do not cross the critical 
line, solutions which start out in this region of the phase plane are thus 
candidates for spatially flat bouncing Universes. Solutions below the critical 
line have ip < and will hence not exhibit a bounce. 

There is a separatrix line between phase plane trajectories which start 
near the origin and which reach ip = 2 (and are thus candidates for a bounc- 
ing Universe) and those trajectories which cross the critical line tp c (H) and 
turn around, i.e. become solutions with ip < 0. The precise location of 
this separatrix is relevant to the initial condition problem of pre-big-bang 
cosmology (see e.g. [22]). The analysis of [20] yields the estimate 

\H\ = -ip 2 withe > 4V2. (29) 



with c S> 4 for the location of the separatrix near the origin of phase space. 

The only region of phase space where singular solutions might occur are 
the critical lines mentioned above. However, one can show [20] that these 
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critical lines not attractors, and hence trajectories peel away from these 
lines as \H\ — > oo and tend to the asymptotic de Sitter region. 

The equations of motion (20) were solved numerically. The resulting 
phase diagram is shown in Figure 2. The absence of singular solutions is 
manifest. In the asymptotic regions \ip\ » 4, all solutions tend to de Sitter 
space. The critical lines are seen to repel contracting solutions towards the 
asymptotic de Sitter regions. The most interesting class of solutions arc 
those corresponding to a cosmological bounce. As predicted, they form a 
set of finite measure among solutions which start out close to the origin of 
the phase plane. There are also solutions which "oscillate" about Minkowski 
space-time. 

In the presence of the dilaton, the phase space becomes three dimen- 
sional (Y>(i), H(t), x(t)) an d therefore more difficult to discuss analytically. 
Note that in the Einstein frame the dilaton corresponds to a homogeneous 
free massless scalar field. Thus, it is already clear from the analysis of [18] 
that the dilaton will not introduce any singularities into the system. In 
fact, from the equation of motion for H (see (20)) it follows that, in the 
region of large \H\, the presence of x W1U n °t change the phase space tra- 
jectories projected onto the (ip, H) plane, a manifestation of the property 
of asymptotic freedom mentioned above. However, the presence of x will 
greatly accelerate the time evolution of ip on the given (ip, H) trajectory. 
This is easy to see for the contracting de Sitter solutions, since in this case 
the x equation of motion (see (20)) leads to exponential growth of x which, 
inserted into the ip equation of motion, demonstrates that at large values 
of | ip\, the x 2 term dominates the evolution of ip. 

Consider now the effect of x on the trajectories in the (ip, H) phase 
plane. The role of the lines ip = ±2 remains unchanged: they correspond 
to maxima of \H\ for any given trajectory. However, the condition for the 
bounce changes. Instead of V(ipb) = 0, it now follows from (20) that the 
condition becomes 

\x 2 + V(i> b ) = 0. (30) 

Hence, \ipb\ is shifted to a larger (and dilaton-dependent) value. Note that 
for very large large values of |%| for which the above equation has no solu- 
tions there will be no bounce. 

The presence of x a ls° changes the critical lines. For a fixed value of ip, 
the condition ip = which determines the critical line occurs at a larger 
value of \H\ than it does in the absence of x> as can easily be seen from 
the ip equation of motion (see (20). In addition, for fixed initial {ip,H) the 
value of ip for contracting solutions is larger with % ^ than with x = 0. 
It thus follows that adding a small value of x increases the range of initial 
conditions in the (ip, H) plane which lead to a bounce. 
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It can be argued that collapsing solutions in the presence of a small \x\ 
quite generically lead to a spatially flat bouncing Universe. Consider initial 
conditions with small but positive x which in the (ip, H) plane lie below the 
separatrix line and which initially "oscillate" about Minkowski space-time. 
These trajectories do not bounce upon reaching H = but start another 
cycle with H < 0. Since H < at all times, %(i) 1S increasing. Eventually, 
x(t) reaches a sufficiently large value such that the trajectory crosses the 
"separatrix" in the (ip, H) plane and evolves past ip = 2 to a successful 
bounce. An example of such a trajectory is shown in Figure 3. The top left 
panel shows the evolution of H as a function of time, the top right depicts 
that of a(t), the bottom left that of x(£) and the bottom right that of ip(t). 

It is now obviously possible to construct initial conditions in which the 
Universe is initially dominated by the dilaton and contracting towards a 
bounce. For example, take the initial conditions for (tp, H, x) to be those of 
the solution in Figure 3 at the end of the cycle preceding the bounce. Thus, 
the nonsingular Universe construction provides a successful implementa- 
tion of the evolution postulated in pre-big-bang cosmology. Note that after 
the bounce, the dilaton rapidly becomes irrelevant to the evolution in the 
expanding phase. Note, in particular, that the dilaton tends to a constant. 
This feature is different from what happens in other attempts to achieve a 
branch change where the dilaton continues to grow in the expanding phase. 

2.5. DISCUSSION 

It has been shown that a class of higher derivative extensions of Einstein 
gravity exists [17, 18] for which all homogeneous and isotropic cosmolog- 
ical solutions are nonsingular. This class of models is very special. Most 
higher derivative theories of gravity have, in fact, much worse singularity 
properties than the Einstein theory. What is special about our class of the- 
ories is that they are obtained using a well motivated Lagrange multiplier 
construction which implements the limiting curvature hypothesis. 

The construction of [17, 18] immediately carries over to dilaton gravity 
[20] and the resulting action solves the graceful exit problem of pre-big- 
bang cosmology, giving automatically a spatially flat bounce interpolating 
between the contracting dilaton-dominated epoch and the expanding FRW 
phase. From the point of view of string theory and pre-big-bang cosmology 
it is more natural to work in the string frame, rather than in the Einstein 
frame (as was done in [20]). However, the nonsingular Universe construction 
can also be applied to dilaton gravity in the string frame, with identical 
results [21]. 

It is to be expected that the nonsingular Universe construction with the 
invariant I2 given by (14) also will regulate the singularities of black holes 
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and of anisotropic homogeneous cosmologies. However, with the exception 
of the case of the two-dimensional black hole [19] the analysis has not been 
completed. 

By construction, all solutions are de Sitter at high curvature. Thus, the 
theories automatically have a period of inflation (driven by the gravity sec- 
tor in analogy to Starobinsky inflation [23]) in the early Universe. A very 
important property of these theories is asymptotic freedom. This means 
that the coupling between matter and gravity goes to zero at high cur- 
vature, and might lead to an automatic suppression mechanism for scalar 
fluctuations. Thus, the nonsingular Universe construction has the potential 
of solving two of the problems of potential-driven inflationary cosmology 
listed in Section 1. 

3. Elements of String Cosmology 

3.1. MOTIVATION 

In the previous section we studied effective actions for the space-time metric 
and the dilaton fields which might arise in the intermediate energy regime 
of string theory. However, it is also of interest to explore the predictions 
of string theory which depend specifically on the "stringy" aspects of the 
theory and which are lost in any field theory limit. It is to a description of 
a few of the string-specific cosmological aspects to which we turn in this 
section. 

3.2. IMPLICATIONS OF TARGET SPACE DUALITY 

Target space duality [24] is a symmetry specific to string theory. As a simple 
example, consider a superstring background in which all spatial dimensions 
are toroidally compactified with equal radii. Let R denote the radius of the 
torus. 

The spectrum of string states is spanned by oscillatory modes which 
have energies independent of R, by momentum modes whose energies E n 
(with integer n) are 

En = £ , (31) 

and by winding modes with energies E' m (m integer) 

E' m = mR. (32) 

Target space duality is a symmetry between two superstring theories, 
one on a background with radius R, the other on a background of radius 
l/R, under which winding and momentum modes are interchanged. 
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Target space duality has interesting consequences for string cosmology 
[4]. Consider a background with adiabatically changing R(t). While R(t) S> 
1, most of the energy in thermal equilibrium resides in the momentum 
modes. The position eigenstates \x > are defined as in quantum field theory 
in terms of the Fourier transform of the momentum eigenstates \p > 

\x >=Y,e ix - p \p> . (33) 
p 

However, for R(t) <C 1, most of the energy flows into winding modes, and 
it takes much less energy to measure the "dual distance" \x > than \x >, 
where 

\x>=J2e i£ - w \w> (34) 

w 

is defined in terms of the winding modes \w >. 

We conclude that target space duality in string theory leads to a min- 
imum physical length in string cosmology. As R(t) decreases below 1, the 
measured length starts to increase again. This could lead to a bouncing or 
oscillating cosmology [4]. 

It is well known that for strings in thermal equilibrium there is a max- 
imal temperature, the Hagedorn temperature [25]. Target space duality 
implies that in thermal equilibrium the temperature in an adiabatically 
varying string background begins to decrease once R(t) falls below 1: 

T (^) = T(R) . (35) 

Thus, the T(R) curve in string cosmology is nonsingular and very different 
from its behavior in standard cosmology. For further discussions of the 
thermodynamics of strings see, e.g., [26] and [27] and references therein. 

3.3. STRINGS AND SPACE-TIME DIMENSIONALITY 

Computations [4] using the microcanonical ensemble show that for all spa- 
tial directions compactified at large total energy E, the entropy S is pro- 
portional to E: 

S = [3 H E, (36) 

with (3h denoting the inverse of the Hagedorn temperature T#. Thus, the 
E{R) curve in string cosmology is very different from the corresponding 
curve in standard cosmology. 

For large R S> 1, most of the energy in a gas of strings in thermal 
equilibrium will flow into momentum modes, and the thermodynamics will 
approach that of an ideal gas of radiation for which 

E(R) ~ I . (37) 
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By duality, for small R 

E{R)~R. (38) 

If, however, for some reason the string gas falls out of equilibrium, the 
E(R) curve will look very different. Starting at R = 1 with a temperature 
approximately equal to Th, a large fraction of the energy will reside in 
winding modes. If these winding modes cannot annihilate, thermal equilib- 
rium will be lost, and the energy in winding modes will increase linearly in 
R, and thus for large R: 

E{R)~R. (39) 

Newtonian intuition indicates that out of equilibrium winding modes 
with an energy relation (38) will prevent the background space from ex- 
panding [4]. The equation of state corresponding to a gas of straight strings 
is 

P = ~JjP ( 4 °) 

where p and p denote pressure and energy density, respectively, and N is the 
number of spatial dimensions. According to standard general relativity, an 
equation of state with negative pressure will lead to more rapid expansion 
of the background. It turns out that the Newtonian intuition is the correct 
one and that general relativity gives the wrong answer [28, 29]. At high 
densities, the specific stringy effects - in particular target space duality - 
become crucial. 

The Einstein action violates duality. In order to restore duality, it is 
necessary to include the dilaton in the effective action for the string back- 
ground. Recall that the action for dilaton gravity can be written as 

S = j d N+1 x^e-^[R + A{V<p) 2 } (41) 

where <fi is the dilaton (rescaled compared to (9)). It is convenient to use 
new fields (p and A defined by 

a(t) = e xt (42) 

and 

ip = 2(f> -NX. (43) 
The action (41) has the duality symmetry 



A — > —A, (p — ► (p . 



(44) 
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The variational equations of motion derived from (41) for a homoge- 
neous and isotropic model are [28, 29] 

2 = e^E + NX 2 (45) 

X-ipX = hsPP (46) 

(p = ^E + NX 2 , (47) 

where P and E are total pressure and energy, respectively. For a wind- 
ing mode-dominated equation of state (and neglecting friction terms) the 
equation of motion for X(t) becomes 

X = ~ef>E(\), (48) 

which corresponds to motion in a confining potential. Hence, winding modes 
prevent the background toroidal dimensions from expanding. 

These considerations may be used to put forward the conjecture [4] that 
string cosmology will single out three as the maximum number of spatial 
dimensions which can be large (R S> 1 in Planck units). The argument 
proceeds as follows. Space can, starting from an initial state with R ~ 1 
in all directions, only expand if thermal equilibrium is maintained, which 
in turn is only possible if the winding modes can annihilate. This can only 
happen in at most three spatial dimensions (in a higher number the prob- 
ability for intersection of the world sheets of two strings is zero). In the 
critical dimension for strings, N = 3, the evolution of a string gas has been 
studied extensively in the context of the cosmic string theory (see e.g., 
Refs. [30], [31] and [32] for recent reviews). The winding modes do, indeed, 
annihilate, leaving behind a string network with about one winding mode 
passing through each Hubble volume. Thus, in string cosmology only three 
spatial dimensions will become large whereas the others will be confined to 
Planck size by winding modes. Crucial aspects of this scenario have been 
confirmed in numerical simulations with cosmic strings in three and four 
spatial dimensions [33]. 

4. A varying speed of light 

One of the most controversial alternatives to inflation is the varying speed 
of light (VSL) theory. This idea was initially proposed by Moffat [6], and 
has recently attracted a lot of attention [5, 34, 35, 36, 37, 38, 39, 40, 41]. 
In the same way inflation is not a model, but a vast class of models, the 
varying speed of light theory may also receive very different realizations. 
Here we present some possible VSL models. 
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4.1. THE INSPIRATION 

The insight leading to VSL concerns the horizon problem. At any given time 
any observer can only see a finite region of the Universe, with comoving 
radius = crj, where rj denotes conformal time, and c the speed of light. 
Since the horizon size increases with time we can now observe many regions 
in our past light cone which are causally disconnected, that is, outside each 
others' horizon (see Fig. 4). The fact that these regions have the same 
properties (eg. Cosmic Microwave background temperatures equal to a few 
parts in 10 ) is puzzling as they have not been in physical contact. This is 
a mystery one may simply relegate to the setting up of initial conditions in 
our Universe. 

It is immediately obvious that the horizon problem may be solved if 
light travelled faster in the Early Universe. Suppose there was a "phase 
transition" at time t c when the speed of light changed from c~ to c + . 
Our past light cone intersects t = t c at a sphere with comoving radius 
r = c + (i]o— rjc), where 770 and r] c are the conformal times now and at t c . This 
is as much of the Universe after the phase transition as we can see today. 
On the other hand the horizon size at t c has comoving radius = c~rj c . If 
c~/c + S> rjo/rjc, then r <C r^, meaning that the whole observable Universe 
today has in fact always been in causal contact (see Fig. 5). Some simple 
manipulations show that this requires 

c~ 1 1 T+ 

logio > 32 - - log 10 z eq + - log 10 y+ (49) 

where z eq is the redshift at matter radiation equality, and T+ and Tp are 
the Universe and the Planck temperatures after the phase transition. If 
T+ ~ Tp this implies light travelling more than 30 orders of magnitude 
faster before the phase transition. It is tempting, for symmetry reasons, 
simply to postulate that c~ = 00. 

This is the initial insight: to replace inflationary superluminal expansion 
by a varying speed of light. The challenge is now to set up a consistent VSL 
theory incorporating this feature, and possibly solving other cosmological 
puzzles, such as the flatness and cosmological constant problems. 

4.2. THE PHYSICAL MEANING OF A VARIABLE SPEED OF LIGHT 

Let us first address the question of the meaning of a varying speed of light. 
Could such a phenomenon be proved or disproved by experiment? Physi- 
cally it does not make sense to talk about constancy or variability of any 
dimensional "constant" . A measurement of a dimensional quantity must al- 
ways represent its ratio to some standard unit. Hence any measurement is 
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necessarily dimensionless. It follows that one can only test experimentally 
the constancy of dimensionless combinations of constants. 

What would we do therefore if we were to observe changing dimension- 
less quantities? Any theory explaining the phenomenon would necessarily 
have to make use of dimensional quantities. It would a priori be a matter of 
choice, prejudice, or convenience to decide which dimensional quantities are 
variable and which are constant (as we mentioned in the illustration above). 
There would be a kind of equivalence, or duality between theories based on 
any two choices as far as dimensionless observations are concerned. How- 
ever, the equations for two theories which are observationally equivalent, 
but which have different dimensional parameters varying, will in general 
not look the same, and again simplicity will end up being an important 
factor in making a choice between theories. 

Consider for instance the recent claim [42] of experimental evidence for 
a time changing fine structure constant a = e 2 /(4ir he). In building a theory 
which explains a variable a we must make a decision. We could postulate 
that electric charge changes in time, or, say, that he must change in time. 
Bekenstein [44] constructs a theory based on the first alternative. He pos- 
tulates a Lorentz invariant action, which does not conserve electric charge. 
Our theory is based on the second choice. We postulate breaking Lorentz 
invariance, a changing he, and consequently non-conservation of energy. 
Any arguments against the experimental meaning of a changing c can also 
be directed at Bekensteins' changing e theory, and such arguments are in 
both cases meaningless. In both cases the choice of a changing dimensional 
"constant" reverts to the postulates of the theory and is not, a priori, an 
experimental issue. 

Given a variable a, and a VSL theory, it must always be possible to 
redefine units so that c and h are constant, and e varies. The two descrip- 
tions should be equivalent with respect of dimensionless quantities. Let us 
assume that measurements of intervals of length dx, time dt, and energy 
dE, are made in the VSL system of units. In this system c oc % oc 1/y/a and 
e = eo- Now define a new system of units such that the same measurements 
in the new system of units lead to results dx, di, and dE such that 



codx 
c\di 
dE 



cdx 
c 2 dt 
dE 



(50) 
(51) 

(52) 



where cq is a constant, to be identified with the fixed speed of light. These 
relations fully specify the new system of units. One may then construct 
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dimensionless ratios in order to identify the constants in the new system: 

(53) 



cdt c dt 
dx dx 



dEdi dEdt 
GdE GdE 



dxc 4 dx c 4 

e 2 e 



dEdx dE dx 
From these we find that in the new system of units 



(54) 

(55) 
(56) 



c = co (57) 

h = h^- = h (58) 

G = G (59) 

e = eo— oc \fa. (60) 

Hence, in the new system, c and h are constants and e varies as y/a. The 
transformation (53) can also be defined by 

dx = dx/e (61) 

dt = dt/e 2 (62) 

dE = dEe 3 (63) 

where e = e/eo is the dielectric constant of the vacuum in the varying-e 
system. 



4.3. THE FOUNDATIONS OF VSL 

In the formulations given to VSL in [5], Lorentz invariance is explictly 
broken. In [45] it has been conjectured that this is a requirement for the 
cosmological problems to be solved in these theories. In its current formu- 
lation the foundations of the theory may be summarized in the following 
postulates. 

Postulate 1. A changing a is to be interpreted as a changing c and % 
in the ratios c oc fi oc a -1 / 2 . The coupling e is constant. 

This postulate merely sets up the theoretical interpretation of the pos- 
sible experimental fact that a changes, in terms of variable dimensional 
quantities. This is a matter of convention and not experiment, as much as 
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a constant He is a matter of convention. With the above choice a system of 
units for mass, length, time, and temperature is unambiguously defined. 

Postulate 2. There is a preferred frame for the laws of physics. This 
preferred frame is normally suggested by the symmetries of the problem, or 
by a criterium such as c = c(t) . 

If c is variable, Lorentz invariance must be broken. Even if one writes 
Lorentz invariant looking expressions these do not transform covariantly. 
In general this boils down to the explicit presence of c in the operator d^. 
Once one admits that Lorentz invariance must be explicitly broken then 
a preferred frame must exist to formulate the laws of physics. These laws 
are not invariant under frame transformation, and one may expect that a 
preferred frame exists where these laws simplify. 

Postulate 3. In the preferred frame one may obtain the laws of physics 
simply by replacing c in the standard (Lorentz invariant) action, wherever 
it occurs, by a field c = c(x^). 

This is the principle of minimal coupling. Because the laws of physics 
cannot be Lorentz invariant it will not hold in every frame. Hence the appli- 
cation of this postulate depends crucially on the previous postulate supply- 
ing us with a favoured frame. This principle may apply in Minkowski space 
time electrodynamics, scalar field theory, etc, in which case the frame in 
which c = c(t) is probably the best choice. The cosmological frame, endowed 
with the cosmic proper time is probably the best choice in a cosmological 
setting. 

Postulate 4 The dynamics of c must be determined by an action princi- 
ple deriving from adding an extra term to the Lagrangian which is a function 
of c only. 

We merely specify that no fields (including the metric) must be present 
in this extra term because we wish minimal coupling to propagate into the 
Einstein's equations. On the other hand it is easy to add to the c Lagrangian 
a potential forcing c to do anything we want. Later we shall explore some 
possibilities. Here we shall simply consider the situation is which c is forced 
to undergo a phase transition in the early Universe. If the c Lagrangian is 
Brans Dicke type, then we have c oc a n as in [37]. We shall explore these 
scenarios further in Section 6. 

4.4. COSMOLOGICAL EQUATIONS 

In a cosmological setting the postulates proposed imply that Friedman 
equations remain valid even when c / 0: 




(64) 
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a 



4irG 




) 



(65) 



a 



3 



where, we recall, pc 2 and p are the energy and pressure densities, K = 
0, ±1 and G the curvature and the gravitational constants, and the dot 
denotes a derivative with respect to proper time. If the Universe is radiation 
dominated, p = pc 2 /3, and we have as usual a oc t 1 / 2 . We have assumed 
that a frame exists where c = c(t), and identified this frame with the 
cosmological frame. 

The assumption that Einstein's equations remain unaffected by decel- 
erating light carries with it an important consequence. Bianchi identities 
apply to curvature, as a geometrical identity. These then imply stress en- 
ergy conservation as an integrability condition for Einstein's equations. If 
c / 0, however, this integrability condition is not stress energy conserva- 
tion. Source terms, proportional to c/c, come about in the conservation 
equations. 

Although this is a general remark we shall be concerned mostly with 
violations of energy conservation in a cosmological setting. Friedman equa- 
tions can be combined into a "conservation equation" with a source term 
in c/c: 



In a flat Universe {K = 0) a changing c does not violate mass conservation. 
Energy, on the other hand, is proportional to c 2 . If, however, K ^ not 
even mass is conserved. 

5. The cosmological problems and VSL 

Even though it is obvious that a varying speed of light theory solves the 
horizon problem, it is far from obvious that other cosmological porblems 
may be solved as well. Here we concentrate on the flatness, cosmological 
constant, and homogeneity problems. 

5.1. THE FLATNESS PUZZLE 

The flatness puzzle can be illustrated as follows. Let p c be the critical 
density of the Universe: 



that is, the mass density corresponding to K = for a given value of a/a. 
Let us define e = 17 — 1 with f2 = pj p c . Then 




(66) 




(67) 




(68) 



22 



ROBERT H. BRANDENBERGER AND JOAO MAGUEIJO 



If p = wpc (with w = 0), using Eqns.(64), (65), and (66) we have: 

P = _ 3 ^(i +U ;) + 2£^ (69) 
p a c 1 + e 

El = _£(2 + (l + e )(l + 3u;)) (70) 
Pc a 

and so 

e = (l + e)e-(l + 3u;) + 2-e (71) 
a c 

In the SBB e grows like a 2 in the radiation era, and like a in the matter 
era, leading to a total growth by 32 orders of magnitude since the Planck 
epoch. The observational fact that e can at most be of order 1 nowadays 
requires that either e = strictly, or an amazing fine tuning must have 
existed in the initial conditions (e < 10~ 32 at t = tp). This is the flatness 
puzzle. 

The e = solution is in fact unstable for any matter field satisfying the 
strong energy condition 1 + 3w > 0. Inflation solves the flatness problem 
with an inflaton field which satisfies 1 + 3w < 0. For such a field e is driven 
towards zero instead of away from it. Thus inflation can solve the flatness 
puzzle. 

As Eqn. 71 shows a decreasing speed of light (c/c < 0) would also 
drive e to 0. If the speed of light changes in a sharp phase transition, 
with \c/c\ 3> a/a, we can neglect the expansion terms in Eqn. 71. Then 
e/e = 2c/c so that e oc c 2 . A short calculation shows that the condition (49) 
also ensures that e<l nowadays, if e « 1 before the transition. 

The instability of the K 7^ Universes while c/c < can be expected 
simply from inspection of the non conservation equation Eq. (66). Indeed 
if p is above its critical value, then K = 1, and Eq. (66) tells us that mass 
is taken out of the Universe. If p < p c , then K = — 1, and then mass is 
produced. Either way the mass density is pushed towards its critical value 
p c . In contrast with the Big Bang model, during a period with c/c < only 
the K = Universe is stable. 

We have assumed in the previous discussion that we are close, but not 
fine-tuned, to flatness before the transition. It is curious to note that this 
need not be the case. Suppose instead that the Universe acquires "natural 
initial conditions" (eg. e ~ 1) well before the phase transition occurs. If 
such Universes are closed they recollapse before the transition. If they are 
open, then they approach e = — 1. This is the Milne Universe, which is our 
case (constant G) may be seen as Minkowski space-time. Such a curvature 
dominated Universe is essentially empty, and a coordinate transformation 
can transform it into Minkowski space-time. Inflation cannot save these 
empty Universes, as can be seen from Eqn. 71. Indeed even if 1 + 3w < 
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the first term will be negligible if e ~ —1. This is not true for VSL: the 
second term will still push an e = — 1 Universe towards e = 0. 

Heuristically this results from the fact that the violations of energy 
conservation responsible for pushing the Universe towards flatness do not 
depend on there being any matter in the Universe. This can be seen from 
inspection of Eqn. (66). 

In this type of scenario it does not matter how far before the transition 
the "initial conditions" are imposed. We end up with a chaotic scenario in 
which Darwinian selection gets rid of all the closed Universes. The open 
Universes become empty and cold. In the winter of these Universes a phase 
transition in c occurs, producing matter, and leaving the Universe very fine 
tuned, indeed as an Einstein deSitter Universe (EDSU). This point seems 
to have been missed in [46]. 



5.2. THE COSMOLOGICAL CONSTANT PROBLEM 

There are two types of cosmological constant problems, and we wish to 
start our discussion by differentiating them. Let us write the action as: 

S = J ^ 4 ^( c4( f 6 ^ Al) + C M + £ A2 ) (72) 

where Cm is the matter fields Lagragian. The term in Ai is a geometrical 
cosmological constant, as first introduced by Einstein. The term in A2 rep- 
resents the vacuum energy density of the quantum fields [47]. Both tend 
to dominate the energy density of the Universe, leading to the so-called 
cosmological constant problem. However they represent two rather differ- 
ent problems. We shall attempt to solve the problem associated with the 
first, not the second, term. Ususally one hopes that the second term will be 
cancelled by an additional couter-term in the Lagrangian. In the rest of this 
paper it is the geometrical cosmological constant that is under scrutiny. 

If the cosmological constant A / then the argument in the previous 
section still applies, with p = p m + pa, where p m is the mass density in 
normal matter, and 

Ac 2 

is the mass density in the cosmological constant. One still predicts Jl m + 
SI A = 1, with Q m = p m /p c and ^a = p\/p c - However now we also have 

a ( Vm\ . 3Kc 2 c 
Pm + 3 a( Pm + ^ =~ PA + ^G^c (74) 
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If A is indeed a constant then from Eq. (73) 

^ = 2°- (75) 

PA C 

If we define eA = pa/ Pm we then find, after some straightforward algebra, 
that 

e A = e A (3*(l + w) + 2 C -±±^) (76) 
\ a c 1 + e / 

Thus, in the SBB model, eA increases like a 4 in the radiation era, like a 3 in 
the matter era, leading to a total growth by 64 orders of magnitude since 
the Planck epoch. Again it is puzzling that eA is observationally known 
to be at most of order 1 nowadays. We have to face another fine tuning 
problem in the SBB model: the cosmological constant problem. 

If c = the solution e\ = is in fact unstable for any w > — 1. Hence 
violating the strong energy condition 1 + 3w > would not solve this 
problem. Even in the limiting case w = —1 the solution ca = is not an 
attractor: e\ would merely remain constant during inflation, then starting 
to grow like a 4 after inflation. Therefore inflation cannot "explain" the 
small value of eA, as it can with e, unless one violates the dominant energy 
condition w > — 1. 

However, as Eqn. (76) shows, a period with c/c <C would drive e\ 
to zero. If the speed of light changes suddenly (\c/c\ 3> a /a) then we can 
neglect terms in a/ a, and so 

£a 2°-J- (77) 



£a(1+£a) cl + e 
which when combined with e/e = 26/ c leads to 

oc -L_ (78) 
l + e A 1 + e V ' 

The exact constraint on the required change in c depends on the initial 
conditions in e and e\. In any case once both e 1 and e\ ps 1 we have 
e\ oc c 2 . Then we can solve the cosmological constant problem in a sudden 
phase transition in which 

logio C — » 64 - \ log 10 z eq + 2 log 10 ^ (79) 

This condition is considerably more restrictive than (49), and means a 
change in c by more than 60 orders of mag nitude, if T+ « Tp . Note that 
once again a period with G/G would not solve the cosmological constant 
problem. 
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As in the case of the flatness problem we do not need to impose "natural 
initial conditions" (ca 1) just before the transition. These could have 
existed any time before the transition, and the argument would still go 
through, albeit with a rather different overall picture for the history of the 
Universe. 

If e\ « 1 well before the transition, then the Universe soon becomes 
dominated by the cosmological constant. We have inflation! The curva- 
ture and matter will be inflated away. We end up in a de-Sitter Universe. 
When the transition is about to occur it finds a flat Universe (e = 0), 
with no matter (p m = 0), and with a cosmological constant. If we rewrite 
Eqn.(76) in terms of e m = p m /pA, for e = and \c/c\ 3> a/a, we have 
£ m = — 2(c/c)(l + e m ). Integrating leads to 1 + e m oc c -2 . We conclude that 
we do not need the presence of any matter in the Universe for a VSL tran- 
sition to convert a cosmological constant dominated Universe into a EDSU 
Universe full of ordinary matter. This can be seen from Eqns. (74)-(75). 
A sharp decline in c will always discharge any vacuum energy density into 
ordinary matter. 

We stress the curious point that in this type of scenario the flatness 
problem is not solved by VSL, but rather by the period of inflation preced- 
ing VSL. 

5.3. THE HOMOGENEITY OF THE UNIVERSE 

Solving the horizon problem by no means guarantees solving the homogene- 
ity problem, that is, the uncanny homogeneity of the currently observed 
Universe across many regions which have apparently been causally discon- 
nected. Although solving the horizon problem is a necessary condition for 
solving the homogeneity problem, in a generic inflationary model solving 
the first causes serious problems in solving the latter. Early causal contact 
between the entire observed Universe allows equilibration processes to ho- 
mogenize the whole observed Universe. It is crucial to the inflation picture 
that before inflation the observable universe in well inside the Jeans length, 
and thus equilibrates toward a homogeneous state. However no such process 
is perfect, and small density fluctuations tend to be left outside the Hubble 
radius, once the Universe resumes its standard Big Bang course. These fluc- 
tuations then grow like a 2 during the radiation era, like a during the matter 
era, usually entailing a very inhomogeneous Universe nowadays. This is a 
common flaw in early inflationary models [48] which requires additional 
fine-tuning to resolve. 

In order to address this problem in VSL a rather complex calculation 
needs to be carried out. It consists in taking the perturbed Einstein's equa- 
tions assuming a changing c, and from them deriving the equations of mo- 
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tions for the matter perturbation. The basic result is that the comoving 
density contrast A and gauge-invariant velocity v are subject to the equa- 
tions: 



A' — ( Sw 1 — ) A = -(1 + w)kv - 2—wU T (80) 

■ a c J a 

a c J \ 1 + w 2k a \ a c J I 



kc 2 w„ ( 2/3 3 (d\ 2 \ 



+- r- kc [—!— + -— - \wU T (81) 

1+w \l + w k 2 c 2 \ a J J v ' 

where k is the wave vector of the fluctuations, and T is the entropy pro- 
duction rate, Tlx the anisotropic stress, and c s the speed of sound given 
by: 

2 _ P' „2^ n 2 1 c' a 



s ft V 31 + wca'J y ' 

In the case of a sudden phase transition Eqn. (80) shows us that Aocc, 
regardless of the chosen equations of state for T and 11^. Hence 

A+ c+ 

A= " ~ (83) 

meaning a suppression of any fluctuations before the phase transition by 
more than a factor of 10 -60 if condition (79) is satisfied. The suppression 
of fluctuations induced by a sudden phase transition in c can be intuitively 
understood in the same fashion as the solution to the flatness problem. 
Mass conservation violation ensures that only a Universe at critical mass 
density is stable, if c/c <C 0. But this process occurs locally, so after the 
phase transition the Universe should be left at critical density locally. Hence 
the suppression of density fluctuations. 

We next need to know what are the initial conditions for A and v. 
Suppose that at some very early time t{ one has c/c = and the whole 
observable Universe nowadays is inside the Jeans length: 770 *C Cir\ij\f?>. 
The latter condition is enforced as a byproduct of solving the horizon prob- 
lem. The whole observable Universe nowadays is then initially in a thermal 
state. What is more each portion of the Universe can be described by the 
canonical ensemble and so the Universe is homogeneous apart from thermal 
fluctuations [49] . These are characterized by the mass fluctuation 



_ 2 _ (5M 2 ) = 4h?i 
(M) 2 Mc 2 



°M = -7T7PT = -T77T ( 84 ) 
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Converted into a power spectrum for A this is a white noise spectrum with 
amplitude 

P A (k) = (\A(k) 2 \) oc ^ (85) 
Pi c i 

What happens to a thermal distribution, its temperature, and its fluctu- 
ations, while c is changing? In thermal equilibrium the distribution function 
of particle energies is the Planck distribution P(E) = \j(e E l khT — 1), where 
T is the temperature. When one integrates over the whole phase space, one 
obtains the bulk energy density pc 2 oc (kbT) 4 / (he) 3 . Let us now consider 
the time when the Universe has already flattened out sufficiently for mass 
to be approximately conserved. To define the situation more completely, we 
make two additional microphysical assumptions. Firstly, let mass be con- 
served also for individual quantum particles, so that their energies scale like 
E oc c 2 . Secondly, we assume particles' wavelengths do not change with c. If 
homogeneity is preserved, indeed the wavelength is an adiabatic invariant, 
fixed by a set of quantum numbers, eg: A = L/n for a particle in a box of 
size L. 

Under the first of these assumptions a Planckian distribution with tem- 
perature T remains Planckian, but T oc c 2 . Under the second assumption, 
we have A = 2irhc/E, and so Ti/c should remain constant. Therefore the 
phase space structure is changed so that, without particle production, one 
still has pc 2 oc (fc{,T) 4 /(/ic) 3 , with Toe c 2 . A black body therefore remains 
a black body, with a temperature T oc c 2 . If we combine this effect with 
expansion, with the aid of Eqn. (66) we have 

T + t(--2-)=0 (86) 
\a cj 

We can then integrate this equation through the epoch when c is changing 
to find the temperature Tj of the initial state. This fully fixes the initial 
conditions for scalar fluctuations, by means of (85). 

In the case of a sudden phase transition we have T + = T~c 2+ /c 2 ~ , and 

so 

2 _ _ Ak b T- 4k b T+ 

" Mc 2 ' ~ Mc 2 + 180 

or 

A-W 2 -^ (88) 

but since A oc c we have 



p^cr^ c 
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Even if T + = Tp = 10 19 Gev these fluctuations would still be negligible 
nowadays. Therefore although the Universe ends up in a thermal state after 
the phase transition, its thermal fluctuations, associated with the canonical 
ensemble, are strongly suppressed. 

6. A possible experimental test of VSL 

VSL, like inflation, may always evade conflict with experiment. In the phase 
transition scenario of Albrecht and Magueijo, changes in c occur only in 
the very Early Universe. They are beyond a chance of contradiction with 
experiment. 

This situation is not entirely satisfactory. A theory which cannot be 
falsified, cannot also be true! In this section we show how it is possible to 
put VS1 to the test. We follow recent work by Barrow and Magueijo [51]. 
In this work, the authors propose a dynamical varying c, that is a theory in 
which the changes in c can be predicted from the matter distribution, rather 
than being just postulated. The model then allows for a joint explanation 
of two recent observations. 

One is the Supernovae Cosmology Project and the High-z Supernova 
Search ([52, 53, 54, 55]). These extended the reach of the Hubble diagram 
to high redshift and provided new evidence that the expansion of the uni- 
verse is accelerating. This may imply that there exists a significant positive 
cosmological constant, A. The other observation [42] was already discussed 
above: recent evidence for a redshift dependence of the fine structure con- 
stant. 

In this section we shall first reinterpret the Hubble diagram in VSL the- 
ories, following which we shall show how a dynamical model for a changing 
c may explain these two observations. 

6.1. THE VSL HUBBLE DIAGRAM 

The Hubble diagram is a plot of luminosity distance against redshift. The 
purpose is to map the expansion factor a(t), where t is the comoving proper 
time. Redshifts provide a measurement of a at the time of emission. If the 
objects under observation are "standard candles" (as Type la supernovae 
are assumed to be), their apparent brightness gives their (luminosity) dis- 
tance, which, if we know c, tells us their age. By looking at progressively 
more distant objects we can therefore map the curve a(t). 

We now examine how this construction is affected by a changing c. 
In [5] we showed that E oc c 2 for photons in free flight. We also showed 
that quantum mechanics remains unaffected by a changing c if fi oc c (in 
the sense that quantum numbers are adiabatic invariants). If H oc c, the 
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Rydberg energy scales like c 2 : 

4 4 

£« = ^2- «c 2 (90) 

Hence all absorption lines, ignoring the fine structure, scale like c 2 . When 
we compare lines from near and far systems we should therefore see no 
effects due to a varying c; the redshift z is still 

l + z e = a /a e (91) 

where o and e label epochs of observation and emission. 

In order to examine luminosity distances, we need to reassess the con- 
cept of standard candles. For simplicity let us first treat them as black 
bodies. Then their temperature scales as T oc c 2 ([5]), their energy density 
scales as p oc T 4 /(hc) 3 oc c 2 , and their emission power as P = p/c oc c, 
implying that standard candles are brighter in the early universe if c < 0. 
However, the power emitted by these candles, in free flight, scales like c; 
each photon's energy scales like c 2 , its speed like c, and therefore its energy 
flux like c. The received flux, as a function of c, therefore scales like: 

P c 2 

P r = -r^T K c ( 92 ) 

where r is the conformal distance to the emitting object, and the subscripts 
r and e label received and emitted. In an expanding universe we therefore 
still have 

where P e Q is the emitting power of standard candles today. Notice that the 
above argument is still valid if the candles are not black bodies; it depends 
only on the scaling properties of emitted and received power. 

We can now set up the Hubble diagram. Consider the Taylor expansion 

a(t) = oo[l + H (t - to) - 2«o#o(* " *o) 2 + •••] (94) 

where Hq = do/ao is the Hubble constant, and qo = —doao/ciQ is the dec- 
celeration parameter. Hence 

z = H {t -*) + (! + q /2)H^(t - t ) 2 + (95) 



or 



to-t=—[z-(l+q /2)z 2 + ...}. (96) 
-no 
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From (93) we find that the luminosity distance di, is 

1 /2 

"^{W -£-*•<>+«.>■ <«7) 
The conformal distance to the emitting object is given by 

r = r*o c (t)dt 

Jt a(t) 

From (94) we have that 

r = c [(t -t) + ^-^H (to-t) 2 + ...} (99) 

where we have assumed that locally c = coa n (that is c = co[l + nHo(t — 
to) + ...]). Substituting (96) we finally have : 

dL = ^-[z + kl-(q + n))z 2 + ...} (100) 

We see that besides the direct effects of VSL upon the expansion rate of the 
universe, it also induces an effective acceleration in the Hubble diagram as 
an "optical illusion" (we are assuming that c decreases in time: n < 0). This 
is easy to understand. We have seen that VSL introduces no intrinsic effects 
in the redshifting spectral line or in the dimming of standard candles with 
distance and expansion. The only effect VSL induces on the construction of 
the Hubble diagram is that for the same redshift (that is, the same distance 
into the past) objects are farther away from us because light travelled faster 
in the past. But an excess luminosity distance, for the same redshift, is 
precisely the hallmark of cosmological acceleration. However, we need to 
consider the other experimental input to our work: the [42] results. By 
measuring the fine structure in absorption systems at redshifts z ~ 0(1) 
we can also map the curve c(t). Since c = co[l + nHo(t — to) + ...] we have 
c = co[l — nz + ...], and so to first order a = ao[l + 2nz + ...]. However, the 
results presented in [42] show that n is at most of order 10~ 5 . This means 
that the direct effects of varying c permitted by the quasar absorption 
system observations are far too small to explain the observed acceleration. 
We need to look at a fully self-consistent generalisation of general relativity 
containing the scope for varying c. 

6.2. THE MODEL 

We start with some general properties of the dynamics of c. Drawing inspi- 
ration from dilaton theories (like Brans-Dicke gravity) we take 

ip = log (c/cq) 



(98) 
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as the dynamical field associated with c. Indeed, powers of c appear in 
all coupling constants, which in turn can be written as e^, where 4> is the 
dilaton. We then endow V with a dynamics similar to the dilaton. The 
left-hand side for the ip equation for a homogeneous field in an expanding 
universe should be V> + 3^ip. This structure ensures that the propagation 
equation for ip is second-order and hyperbolic. However, since VSL breaks 
Lorentz invariance one need not choose (as in Brans-Dicke) the source term 
to be p — 3p, where p and p are the energy density and pressure of matter 
respectively. A possibility is 

^ + 3^ = 4ttGlo P (101) 

We could have added pressure to the source and terms in the curvature. 
An interesting alternative is 

$ + 3-V> = AnGuip - p/c 2 ) - (102) 
a cr 

These dynamical equations, when combined with Friedmann equations, lead 
to a "Machian" scenario, in which c = Coa n . If we consider the dynamical 
equation (102), then with n = uj, regardless of the equation of state, or 
amount of curvature domination. This is the scenario initially envisaged by 
Barrow [34]. We now see that the coupling constant uj, describing the ability 
of matter to drive changes in the speed of light, is also what determines the 
exponent n in the Machian solution. The flatness problem is solved if: 

uj < (1 - 3 7 /2) (103) 

Effectively this means uj < — 1, since the flatness problem should be solved 
in the radiation epoch. This is a self regulating scenario, in which curvature 
drives changes in the speed of light, which in turn suppress curvature. If 
we take Eqn.(lOl) then if the curvature term is subdominant one has 

ra = w/(2-7) (104) 

that is n = 3/2uj in the radiation epoch, and n = uj in the matter epoch. 
The flatness problem is solved if 

oo < (2- 7 )(l-3 7 /2) (105) 

Effectively this means uj < —2/3, since the flatness problem should be 
solved in the radiation epoch. Assuming that curvature has already been 
dominated is natural in these scenarios, as there is no Planck time. Hence 
there is no natural time when we may give matter and curvature initial 
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values of the same order of magnitude. It is sufficient for flatness to be an 
attractor. 

In order to explain the Supernovae results, however, the best dynamical 
equation is 

(106) 



a 



4ttGuj^, 



where p is the total pressure of the matter fields and uj is a coupling constant 
(distinct from the Brans Dicke coupling constant). The full self-consistent 
system of equations in a matter-plus-radiation universe containing a cos- 
mological constant stress 



PA 



is therefore 



4> + 3-ip 
a 

A O 

p 7 + 4-p 7 
a 

PA 

CL 

Pm + 3 — Pm 

a 



-)' 

a) 



Ac 2 
4vrG^, 

-2p A i;, 

tPAip, 

o, 

8ttG 



(Pm + Pj + PA), 



(107) 

(108) 
(109) 
(110) 

(111) 



where subscripts 7 and m denote radiation and matter respectively. We 
have assumed that the sink term (109) is reflected in a source term in 
(108) (and not in (110)). This is due to the fact that this term is only 
significant very early on, when even massive particles behave like radiation. 
We have ignored curvature terms because in the quasi-lambda dominated 
solutions we are about to explore we know that these are smaller than p\, 
([37]). Here, in complete contrast to Brans-Dicke theory, the field ip is only 
driven by radiation pressure in the dust-dominated era. In other words, 
only conformally invariant forms of matter couple to the field ip. 



6.3. SOLUTIONS TO THE MODEL 

In a radiation-dominated universe the behaviour of this system changes at 
the critical value u = —4. For uj < —4 we reach a flat A = attractor as 
t — > 00. For — 4 < uj < we have attractors for which pa and p 7 maintain 
a constant ratio (see [37]). In Fig. 6 we plot a numerical solution to this 
system, with uj = —4.4 (a 10% tuning below the critical value u = —4) and 
n = —2.2 during the radiation epoch. As expected from [37], this forces Oa 
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to drop to zero, while the expansion factor acquires a radiation-dominated 
form, with a oc t 1 / 2 . By the time the matter-dominated epoch is reached, Qa 
is of order 10 -12 . During the matter epoch, the source term for ip disappears 
in eq. (106), n starts to approach zero, Qa starts to increase, and the 
expansion factor takes on the a oc t 2//3 dependence of a matter-dominated 
universe. A few expansion times into the matter epoch, S7a becomes of 
order 1 and the universe begins accelerating. By the time this happens n 
is of order 10 -5 , in agreement with the expectations of [42]. This type of 
behaviour can be achieved generically, for different initial conditions, with 
a tuning of uj that never needs to be finer than a few percent. 

We can provide an approximate argument explaining why this theory 
should display this type of behaviour and why we need so little fine tuning 
of uj to explain the supernovae experiments. If we neglect changes in c after 
matter-radiation equality, t eq , we are going to require 

P^eq) ^ ^-3 _ 1n -l2 



P(teq) 



z-< ~ lO" 12 . (112) 



Let c = coa n( ^\ with n = —2 — 5, and n = uj/2 during the radiation epoch. 
We can integrate the conservation equations to give 

P_ = A !L_ M 13) 

p A a 4 pA n + 2' v ' 

with A constant, from which it follows that 



P_ = 2 

PA 5 



2p\iJ \a 



(114) 



We see that assymptotically p/p\ grows to infinity, if 5 > (the flat A = 
attractor of [37]). However the growth is very slow even if 5 is not very small. 
Our theory displays very long transients, and a very slow convergence to its 
attractor, a property similar to quintessence models ([50]). It is therefore 
possible to achieve pa/p ~ 10 -12 at the end of the radiation epoch, with S 
chosen to be of order 0.1. 

Now, why is the change in c of the right order of magnitude to explain 
the results of [42]? With a solution of the form c = cqo! 1 ^ we find that 

n{t) « - ^ (115) 

V ' 3(p m + 2p A ) V ; 

With uj ~ —4 we therefore have 

4 2.3 x 10~ 5 . . 

" ( '° ) a -swaTsw (U6) 
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of the right order of magnitude. The order of magnitude of the index n ~ 
1CT 5 , observed by [42], is therefore fixed by the ratio of the radiation and 
the matter energy densities today. 

In summary in this theory gravitational effect of the pressure drives 
changes in c, and these convert the energy density in A into radiation. Thus 
A is prevented from dominating the universe during the radiation epoch. 
As the universe cools down, massive particles eventually become the source 
of pressureless matter and create a matter-dominated epoch. In the matter 
epoch the variation in c comes to a halt, with residual effects at z « 1 — 5 
at the level observed by Webb et al. As the c variation is switched off, the 
A stress resurfaces, and dominates the universe for a few expansion times 
in the matter-dominated era, in agreement with the supernovae results. 

7. A bit of wishful thinking 

From all that was said above it should be clear that these ideas require and 
merit further work. Here we summarise what one may expect in the near 
future. 

Obtaining a nonsingular cosmological model remains one of the out- 
standing goals of string cosmology. However, it will be crucial to go beyond 
the effective field theory models discussed in Section 2, and of the per- 
turbative string framework used in Section 3. The recent developments in 
string/M theory may shed a new light on the singularity issue. The new 
fundamental degrees of freedom in string theory associated with d-branes 
will inevitably play a crucial role in cosmology; for the singularity problem, 
the dynamics of expansion and compactification, and for many other issues. 
We should expect major progress in this area in the coming years. 

While VSL scenarios remain an interesting idea in cosmology, it is clear 
that a lot of work on its foundations needs to be done. The model requires 
breaking Lorentz symmetry and general covariance. Clearly this is far from 
desirable. Moffat's early formulation [6], relies on spontaneous, rather than 
explicit, breaking of Lorentz symmetry. The idea is that the vacuum it- 
self is not Lorentz invariant, even though the full theory is. More recently 
Moffat and Clayton [38, 39] and Drummond [41] proposed a bimetric for- 
mulation for VSL. In these models the light-cone of gravity and matter 
are different. These models are on firmer foundations than the so-called 
Albrecht-Magueijo model. Unfortunately their application to cosmology is 
rather cumbersome. In some recent work, it has been suggested that VSL 
might be realized in the brane-world. If the Universe is a 4 dimensional 
brane in the process of being absorbed by a black hole, perhaps the pro- 
jected theory of gravity might resemble the Albrecht and Magueijo theory 
[56, 57]. Notice that in these theories one does not need to postulate a 
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dynamics for the field ip = log(c/co). The dynamics of c follows from the 
kinematics of the brane moving in the space of the black hole. 

An orthogonal approach has been pursued by Manida [58] and Stepanov 
[59]. These authors have noted that breaking Lorentz invariance does not 
imply breaking the principle of relativity; merely that the constancy of the 
speed of light must be abandoned. The symmetry transformation replacing 
Lorentz symmetry in such circumstances is the so-called projective Lorentz 
transformation (or Fock-Lorentz transformation). The ensuing theory al- 
lows for a varying speed of light without favouring one inertial frame over 
another. So far this approach has not incorporated gravity; a necessary 
ingredient before a realistic application to cosmology may be discussed. 

The purpose of [5] was to show how far reaching VSL can be. The 
challenge is now to set up a fundamental theory is which the model might 
be realized. To our mind, the challenge is still open. 
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H 



Figure 1. Penrose diagrams for collapsing Universe (left) and black hole (right) in 
Einstein's theory (top) and in the nonsingular Universe (bottom). C, E, DS and H stand 
for contracting phase, expanding phase, de Sitter phase and horizon, respectively, and 
wavy lines indicate singularities. 
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Figure 2. The phase portrait for solutions of the equations of motion, with the dilaton 
kinetic energy (x) set to zero is shown. H is plotted on the vertical axis, while ip is plotted 
on the horizontal. 
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Figure 3. A specific bouncing solution with a non-trivial contribution from the dilaton 
is plotted. The left top panel shows the evolution of H as a function of time, the right 
top panel depicts that of a(t), on the bottom left the evolution of x(t) is shown and on 
the bottom right that of ip(t). 



TIME 




Figure 4- Conformal diagram (light at 45°) showing the horizon structure in the SBB 
model. Our past light cone contains regions outside each others' horizon. 
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TIME 




Figure 5. Diagram showing the horizon structure in a SBB model in which at time t c 
the speed of light changed from c~ to c + <JC c~ . Light travels at 45° after t c but it travels 
at a much smaller angle with the space axis before t c . Hence it is possible for the horizon 
at t c to be much larger than the portion of the Universe at t c intersecting our past light 
cone. All regions in our past have then always been in causal contact. 
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10 20 28 30 

log(a) log(a) 

Figure 6. Evolution of t x h (where h — ^), n = jj^, and 51a in log or linear plots 
as appropriate. The panels on the left (right) describe the radiation (matter) dominated 
epoch. We have taken ui = —4.4 (a modest 10% tuning over the critical value to = —4). In 
the radiation epoch n = —2.2, 57a slowly drops to zero, and the expansion factor has the 
usual dependence a oc t 1 ^ 2 . As the Universe enters the matter epoch n starts dropping 
towards zero, a oc t 2 ^ :i , and then Qa starts to increase. Eventually n is of order 1CP J and 
Q,a is of order 1. This type of behaviour occurs for a large, non finely tuned, region of 
couplings u) and initial conditions. 



